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SUMMARY 


The  objectives  of  the  project  "Biological -Mathematical  Modeling  of 
Chronic  Toxicity"  were  to  study  the  factors  affecting  the  fate  of  inhaled 
vapors  in  the  body  and  to  develop  a  mathematical  model  to  describe  their 
effects  on  uptake,  distribution  and  elimination. 

The  main  accomplishments  of  the  research  are:  1)  A  general  pharmacoki¬ 
netic  model  for  inhalation  administration,  based  on  physiological  parameters 


of  the  exposed  subject  and  on  physical  and  chemical  properties  of  the  inhaled 


substance,  was  defined.  Mathematical  theories  and  existing  programs  used  in 
electric  engineering  for  solution  of  electric  networks  were  suggested  for  solu¬ 
tion  of  differential  equations  and  specific  programs  were  prepared.  2)  Infor¬ 
mation  on  interspecies  differences  was  extracted  from  the  literature,*  partition 
coefficients  of  eight  organic  solvents  were  determined  for  eight  tissues  of  three 
species  (man,  monkey,  rat);  steady  state  clearances  calculated  from  pulmonary 
uptake,and  intrinsic  clearance  determined  from  vapor  distribution  in  rat  were 
used  in  the  model  as  elimination  rate  constants.  3)  The  models  were  used  to 
evaluate  the  effects  of  the  following  parameters  on  uptake,  distribution  and 
elimination  of  Inhaled  vapors:  solubility,  metabolism,  body  build,  inter-species 
differences,  physical  exertion,  exposure  duration,  exposure  repetition,  and 
short-  term  excursion  limit.  4)  Nonlinear  dependence  of  pulmonary  uptake  on  ex¬ 


posure  concentration  was  observed  in  the  animal  model.  5)  Reduced  pulmonary 


uptake  and  quantitative  and  qualitative  changes  in  elimination  were  observed  if 
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The  objectives  of  the  project  "Biological -Mathematical  Modeling  of 
Chronic  Toxicity"  were  to  study  the  factors  affecting  the  fate  of  inhaled 
vapors  in  the  body  and  to  develop  a  mathematical  model  to  describe  their 
effects  on  uptake,  distribution  and  elimination.  The  ultimate  goal  was  to 
be  able  to  employ  the  mathematical  model  to  predict  optimum  exposure  conditions 
for  biological  testing  of  chronic  toxicity. 

In  our  attempt  to  explore  mathematical  modeling  in  inhalation  toxicology, 
our  research  took  two  courses:  1)  Development  of  an  animal  model  to  generate 
information  for  the  simulation  model  and  to  corroborate  experimentally  the  pre¬ 
dictions  made  by  the  simulation  model.  2)  Development  of  a  general  simulation 
model,  based  on  physiological  parameters,  for  quantitative  description  of  uptake, 
distribution  and  elimination  of  inhaled  vapors  or  gases.  The  general  simulation 
model  was  refined  in  order  to  make  possible  studies  of  various  factors  which  may 
affect  the  fate  of  inhaled  vapors  in  the  body.  The  general  simulation  model  was 
also  simplified  in  order  to  accomodate  the  program  on  a  small  programmable  pocket 
calculator.  The  simplified  model  uses  relatively  inexpensive  equipment,  and  the 
user  needs  no  mathematical  skill.  The  specific  simplified  programs  are  designed 
to  be  used  by  practicing  toxicologists  and  hygienists.  The  program  for  evalua¬ 
tion  of  short-term  excursion  limits  is  an  example. 

The  main  accomplishments  of  the  research  are: 

1)  A  general  pharmacokinetic  model  for  inhalation  administration  was  defined. 
This  model  is  based  on  physiological  parameters  of  the  exposed  subject  and  on 
physical  and  chemical  properties  of  the  Inhaled  substance.  The  concept  of  compart¬ 
ment  model  was  applied. 
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,  2)  Information  on  inter-species  differences  in  pulmonary  ventilation,  tissue 

perfusion,  and  body  build  was  extracted  from  the  literature  to  be  available  for 
substituting  in  the  model. 

3)  Partition  coefficients  of  eight  organic  solvents  were  determined  for  eight 
tissues  of  three  species  (man,  monkey,  rat).  Species  differences  and  effect  of 
meal  Intake  on  partition  coefficients  were  under  scrutiny. 

4)  Steady  state  clearances  calculated  for  man,  dog  and  monkey,  from  pulmonary 
uptake  rate, were  used  to  define  the  over-all  elimination  rate  constants  in  the 
model . 

5)  Intrinsic  clearances  determined  from  vapor  distribution  in  rat  at  steady 
state  were  used  in  the  model  as  elimination  rate  constants  of  individual  elimination 
pathways . 

6)  Mathematical  theories  and  existing  programs  used  in  electric  engineering 
for  solution  of  electric  networks  composed  from  capacitances  and  conductances  are 
suggested  for  solution  of  differential  equations  describing  uptake,  distribution 
and  elimination  of  inhaled  vapors. 

7)  Specific  programs  were  prepared,  the  most  important  being  for:  a)  A  10- 
compartment  model  in  Fortran  IV,  to  be  used  in  large  time-shared  computer  Uni vac 
1100/20.  b)  A  5-compartment  model  in  basic,  to  be  used  on  minicomputer  Apples  II 
Plus,  c)  A  2-compartment  model  to  be  used  in  programmable  pocket  calculator  TI-59. 

8)  The  models  were  used  to  evaluate  the  effects  of  the  following  parameters  on 
uptake,  distribution  and  elimination  of  inhaled  vapors:  a)  solubility,  b)  metabolism, 
c)  body  build,  d)  inter-species  differences,  e)  physical  exertion,  f)  exposure 
duration,  g)  exposure  repetition,  h)short-term  excursion  limit. 

9)  Nonlinear  dependence  of  pulmonary  uptake  on  exposure  concentration  was  ob¬ 
served  In  the  animal  model.  Plateau  kinetics,  applicable  to  capacity-limited  pro¬ 
cesses,  Is  suggested  for  description  of  elimination  processes  in  the  model. 
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10)  Reduced  pulmonary  upake  and  quantitative  and  qualitative  changes  In 
elimination  were  observed  if  two  vapors  were  inhaled  simultaneously. 

Reports  and  Publications; 

(Under  Professional  name  of  P.I.  -  FISEROVA-BERGEROVA) 

Methods  used  in  this  project  and  obtained  results  were  reported  in  detail  in 
interim  scientific  reports  submitted  to  AFOSR  in  March  of  1979  and  1980,  and  pro¬ 
gress  reports  submitted  in  October  of  1978,  1979,  and  1980.  The  final  interim 
scientific  report  follows.  Reprints  or  xerox  copies  of  publications  so  far  re¬ 
sulting  from  this  project  are  enclosed. 

The  theories  and  solutions  of  the  model  and  supporting  data  obtained  during 
this  project  are  being  prepared  for  publication  in  "Modeling  of  the  Uptake  Metabolism 
and  Elimination  of  Some  Vapors  and  Gases",  edited  by  V.  Thomas.  This  monograph  will 
be  published  by  CRC  Press  in  1982.  Manuscripts  of  four  chapters  dealing  with  data 
and  model  developed  and  verified  under  the  AFOSR  contract  will  be  submitted  to 
AFOSR  in  Fall  1981. 

Detailed  Scientific  Report  for  period:  October  1980-May  1981 

The  main  effort  was  concentrated  on  finishing  the  general  program  for  mini¬ 
computer  Apple  II  Plus,  and  on  preparing  manuscripts  for  the  monograph. 

The  program  solves  the  general  multi  compartment  model  for  simulation  of  inha¬ 
lation  administration.  First  order  kinetics  is  assumed.  The  listing  of  the  program 
is  attached. 

The  general  model  denoted  as  LINEAR  has  the  following  options: 

1)  The  model  can  simulate  exposure  of  any  subject  with  defined  physiologic 
parameters  to  any  compound  for  which  partition  coefficients  and  steady  state  clear¬ 
ance  (or  intrinsic  clearance)  are  known.  From  the  mathematical  point  of  view,  clear- 
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ance  can  be  located  in  any  tissue  or  compartment.  The  compartment  can  include 
one  or  several  tissues  with  similar  parameters.  The  compartments  are  numbered 
and  the  numbers  are  used  as  indexes,  (i),  of  the  elements  in  the  model.  The 
scheme  of  the  general  model  is  in  Figure  1.  The  numbering  of  the  elements  (denoted 
by  index  i)  is  important.  The  last  elements  (i=n)  always  relate  to  respiration. 
(G^  =  alveolar  ventilation;  C^=  FRC  +  2/3  +  lung  tissue  volume  multiplied 

by  lung-air  partition  coefficient  +  volume  of  arterial  blood  multiplied  by  blood- 
air  partition  coefficient;  V^  =  concentration  in  alveolar  air  =  concentration  in 
arterial  blood  divided  by  blood-air  partition  coefficients).  It  is  convenient  to 
start  numbering  with  those  organs  or  compartments  in  which  elimination  takes  place, 
and  to  assign  larger  numbers  to  compartments  without  clearance.  The  values  of  the 
elements  are  defined  as  follows:  =  tissue  (or  compartment)  volume  multiplied 

by  tissue-air  partition  coefficient.  G^s  perfusion  multiplied  by  blood-air 
partition  coefficient.  Gx^=  clearance  multiplied  by  blood-air  partition  coeffi¬ 
cient.  If  there  is  no  clearance  in  the  compartment,  Gx(-j)=  AH  partition  coeffi 
cients  relate  to  37°C.  Table  1  is  designed  for  convenient  preparation  of  input 
data  (example  in  table  3). 

FIGURE  1 


General  n-Compartment  Simulation  “odel 


The  element  GX  is  removed  from  those  compartments  in  which  no  clearance 
takes  place. 


TABLE  1 
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TABLE  2 

Comments  to  output  commands 


STATEMENT 

EXPLANATION 

DENOTATION 

1270 

Prints  times  from  the  simulation  start  in  minutes 
(Defined  in  statement  62). 

TIME  = 

1280-1289 

CHOICE:  Concentration  in  tissues  divided  by 
tissue-air  partition  coefficient  Cx,.»3.  Index 
(i)  according  to  compartment  number.  ' 

1280 

Example  in  listing 

X(l) 

X(2) 

1281 

X(3) 

X(4) 

• 

add  additional  compartment  if  .required 

1290-1299 

CHOICE:  Elimination  rates  CCUR/.i].  Index  (i) 
indicates  compartment  number  except  with  the  last 
one.  The  last  CUR  always  indicates  the  uptake 
rate  or  exhalation  rate. 

1290 

Example  in  listing:  CUR/~\  is  either  uptake  rate(+] 
or  exhalation  rate  (-).  ' 

CUR(i) 

CUR(2) 

• 

add  additional  clearances  if  required 

1300-1309 

CHOICE:  Total  amounts  metabolized  or  excreted 
CCYNT^ . j].  Index  indicates  compartment  number 

1300 

Example  in  listing 

CYNT(d 

add  additional  clearances  if  required 

1310 

Total  pulmonary  uptake;  Total  amount  exhaled. 

The  indexes  do  not  relate  to  the  model.  They 
change  automatically  according  to  number  of  excre¬ 
tory  pathways  (defined  in  statement  50). 

Example  in  listing:  CYNT,~v=  uptake(+), 

CYNT(4)»  exhaled(-) 

CYNT(3) 

cynt(4) 
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2)  The  model  accomodates  any  number  of  changes  in  exposure  concentration, 
thus  enabling  simulation  of  any  kind  of  exposure.  In  the  input,  "number  of  time 
section"  (statement  64)  controls  the  number  of  changes  in  exposure  concentrations. 
Information  on  exposure  concentration  and  duration  of  exposure  is  presented  in  con¬ 
secutive  order  in  statements  66-99  (table  1). 

3)  The  output  offers  the  following  information  (tables  2  and  4):  a)  Asympto¬ 
tic  values.  Asymptotic  values,  multiplied  by  appropriate  tissue-air  partition 
coefficient,  represent  tissue  concentrations  reached  at  steady  state  for  exposure 
concentration  equal  to  1.  b)  Time  constants  and  rate  constants,  which  are  hybrid 
constants  related  to  all  distribution  and  elimination  processes.  These  rate  const¬ 
ants  equal  exponent  constants  in  Laplace  transform,  c)  For  time  intervals,  schedul¬ 
ed  in  minutes  (controlled  by  STEP  in  statement  62),  information  is  displayed  in  the 
following  order:  First,  concentration  in  all  compartments  (for  i  =  1  to  n),  divided 
by  partition  coefficients .  Second,  excretion  rates,  or  metabolic  rates,  related 
to  clearance  Gx^.j(for  i  =  1  to  n-1),  followed  by  pulmonary  uptake  rate  (during 
saturation,  positive  sign)  or  exhalation  rate  (during  desaturation,  negative  sign). 
Last,  the  total  amounts  excreted  or  metabolized  via  clearance  denoted  G^.j,  follow¬ 
ed  by  total  pulmonary  uptake  and  total  amount  exhaled. 

4)  The  output  data  is  displayed  on  the  screen.  On  command,  the  output  can 
also  be  printed  or  filed  on  the  disk.  The  data  file  is  arranged  so  that  it  can  be 
used  as  input  in  the  Appleplot  program  for  graphic  display. 

Arrangement  of  input  data  and  output  data  is  shown  in  a  simple  example.  Four 
compartment  model  should  be  used  to  simulate  eight-hour  exposure  of  a  resting  man 
to  typical  organic  solvent  (1  mg/£),  followed  by  16  hour  exposure  to  zero  concentra¬ 
tion.  The  data  for  each  eight  hours  should  be  printed.  The  simulation  model  is  pic¬ 
tured  in  Figure  2.  The  input  data  is  prepared  in  table  3  and  the  printout  with  comment 
ary  is  in  table  4.  The  input  data  is  also  part  of  listing  of  the  program  (framed  state 
ments  in  listing). 
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FIGURE  2 


Simulation  Model  foc_8iHQur_ Exposure  to  Orqanic  Solvents 


tissues 


Data  refers  to  tables  3-4.  Compartment  numbers  are  circled. 


TABLE  3 
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Preparing  Ott*  for  LINEAR 


Simulation  of  8-Hour  Exposure  to  Organic  Solvent* 


STATEMENT 

READ!  IK 

50 

No.  of  coavartaents,  No.  of  clearances 

4 

/ 

52-57 

No.  compartment  :perfus1on)At1S/,l1r,  volume 

FXbl/»1r 

VXbl/e1r 

52 

l:  4i/x/0  .  t  j/x20 

44 

7Z$ 

53 

2.838/0;M^X2O 

2S.3 

U8 

54 

3:  description  /  J2  X  /O  ')  J2. 5 X  400 

4.2 

5000 

55 

4:  description  __ 

— 

— 

58 

n;  alveolar  ventilation,  V,  of  lung* 

/>V*  4 

6 

/£ 

60 

Clearance  In  conp.  ©2,3,4 

0. 847  x  JO 

$.47 

0 

0 

0 

62 

Printer  (Indicate  In  minutes) 

480 

64 

How  nany  exposure  concentrations? 

2 

Indicate  exposur*  conctntntlon  and  duration  tine  In  sequence: 

66 

1:  Exposure  cone.,  duration  (minutes) 

/ 

480 

67 

2: 

0 

no 

68 

3: 

68 

4: 

70 

5: 

90 

n 

f  •  FRC  ♦  2/3  »t1^  .♦  volume  of  lung  tissue x  *  volume  of  »rt.  blood  x  Vl/»1r 


♦Data  are  In  the  listing  of  program  (example). 
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TABLE  4 


Print  out  for  8-hour  exposure  to  1  mg /£  of  an  organic 
sol  vent. 

(Input  data  in  program  listing) 


RSVMPTOTTC  VPL'JES  CCR  UNIT  EXPOSURE 
concentration 


(t)  .334  Vessel  rich  tissue  (elimination! 

g)  .45732  muscle,  skin 

.45792  fat,  marrow 

(|!)  .45792  alveolar  air  or  arterial  blood 


) 


To  obtain  steady  state  concentration  in  tissue  (mg/£),  the  assymptotic  values 
must  be  multiplied  by  appropriate  tissue-air  partition  coefficients. 


STEP*  480 

*  OF  TIME  SECTIONS* 

CONCENTR.  mg/£  TIME  IN  MIN 
1  430 

0  960 


EXPONENT  TIME  CONST,  minute 

-1 . 16583694E-03  S57.752S84 
-.0187986412  53.1953343 

-.246105268  4.06330188 

-7.64555918  .130794881 


concentrations* in 

metabolic  rate 
total  metabolized 
total  uptake 


TIME=480  min 

.274469211 

.186312692 

2.32475422 

911.873557 

2105.08547 


concentrations* in 

metabolic  rate 
total  metabolized 
total  uptake 


-  TIME*9S0  min 
©  .0469573777 
Q)  .116414466 
©  .397728989 

Q  1220.2008 

®>  2105.08547 


TIME* 1448 

.0268251347 

.0665243045 

.227209315 

1366.36937 

2105.08547 


SATURATION 


after  start  of  simulation 
(»  .323863796 
®  .327670372 
(P  4. 03397777  uptake  rate 


®  0 


v 


total  exhaled 


after  start  of  simulation 

§.0574781799 
.0558367587 

-.335020552  exhalation  rate 
<P  -256.65618  total  exhaled 


.0318980546 

-.191388328  exhalation  rate 
-379.78012 


DESATURATION 


*To  obtain  concentrations  In  tissues  (mg/£)  the  values  must  be  multiplied 
by  appropriate  tissue-air  partition  coefficients. 

The  circled  numbers  match  compartment  numbers  in  Table  3  and  Figure  2. 
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GENERAL  DISCUSSION 

Toxic  effect  is  related  to  the  concentration  of  the  substance  in  the  target 
organ  and  to  the  interaction  of  the  substance  with  receptors.  The  concentration 
in  the  target  organ  is  different  from  exposure  concentration  and  changes  with  ex¬ 
posure  duration  in  a  predictable  way.  The  interaction  with  receptors  is  less  pre¬ 
dictable.  The  receptor  can  be  any  constituent  in  the  body,  usually  macromolecules 
including  enzymes  which  intermediate  the  metabolism  of  the  inhaled  substance.  The 
simulation  by  our  mathematical  model  predicts  quantitatively  how  the  volatile  subs¬ 
tance  gets  from  the  environment  to  the  tissues,  but  provides  no  information  on  inter¬ 
action  with  receptor  on  biotransformation,  or  on  toxic  action  mechanism. 

Modern  computer  technology  and  mathematical  theories  enable  solution  of  multi¬ 
compartment  models  with  simultaneous  elimination  processes.  Pharmacokineti cists 
proposed  numerous  models  for  dosing  and  elimination  of  drugs  administered  intravenous¬ 
ly  or  orally.  Anesthesiologists  offered  a  model  for  uptake  of  inert  Inhalation  anes¬ 
thetics  administered  at  constant  alveolar  concentration.  We  proposed  the  first 
multi  compartment  model  for  simulation  of  exposure  to  metabolized  (or  excreted) 
vapors  and  gases  (Fiserova-Bergerova ,V. ,Vlach,  J. ,Singhal ,K. .Simulation  and  pre¬ 
diction  of  uptake,  distribution,  and  exhalation  of  organic  solvents,  Br.J.Ind.Med. , 

31,  45,  1974).  Since  then,  we  refined  the  model  and  searched  for  experimental  de¬ 
signs  suitable  for  obtaining  input  data  for  the  model.  Obtaining  input  data  is  now¬ 
adays  a  more  difficult  task  than  mathematical  solution  of  the  model.  The  following 
information  and  methods  were  missing,  or  were  inadequate,  when  we  prepared  the 
general  simulation  model:  Methodology  for  determination  of  individual  elimination 
rate  constants  In  vivo  is  not  fully  developed-,  diffusion  constants  are  known  for  a 
very  limited  number  of  compounds;  information  on  partition  coefficients  Is  scattered; 
a  non-invaslve  method  for  determination  of  cardiac  output  and  Its  distribution  is 
not  comnonly  available;  the  significance  of  extrahepatic  metabolism  has  only  recent- 


ly  been  appreciated;  no  attempt  has  been  made  to  quantitatively  evaluate  adsorp¬ 
tion  and  desorbtion  of  vapors  in  respiratory  airways.  We  were  able  to  fill  some 
of  these  gaps,  such  as  developing  animal  models  for  determination  of  clearance,  or 
obtaining  data  on  clearance  and  partition  coefficients  of  some  substance.  But  for 
information  on  other  factors,  such  as  diffusion  and  perfusion,  we  relied  only  on 
the  literature.  Some  new  factors  affecting  uptake  and  distribution  of  inhaled 
vapors  we  recognized  during  the  process  of  model  development.  Because  of  lack  of 
time  and  unavailability  of  data,  we  were  not  able  to  include  these  factors  in 
the  model.  Our  observation  concerns  three  particulars: 

1)  Nonlinear  elimination  processes.  The  nonlinear  dependence  of  pulmonary 
uptake  on  exposure  concentration  was  observed  and  was  attributed  to  capacity 
limited  metabolism  or  excretion.  Methods  for  measuring  some  parameters  describing 
nonlinearity  were  developed,  and  a  limited  amount  of  information  was  obtained  in 
our  laboratory  and  elsewhere  (Toxicology  Research  Laboratory  of  Dow  Chemical  Co.;Toxic 
Hazards  Research  Unit  Department  of  the  Air  Force).  Mathematical  theories  suitable 
for  solving  the  nonlinearity  of  elimination  processes  were  scrutinized. 

2)  Additional  factors  affecting  pulmonary  gas  exchange.  We  observed  that 
the  values  of  pulmonary  uptake  rate  obtained  from  the  difference  between  concentra¬ 
tions  in  ambient  air  and  mixed  exhaled  air  multiplied  by  minute  ventilation  were  in 
some  instances  smaller  than  the  values  obtained  from  the  difference  between  concen¬ 
trations  in  ambient  air  and  end  exhaled  air  multiplied  by  alveolar  ventilation.  We 
found  three  possible  explanations:  metabolism  of  vapor  in  lung  tissue;  adsorbtion 
and  desorbtion  of  vapors  in  respiratory  airways;and  slow  diffusion  accross  alveolar 
membranes.  The  first  possibility  is  supported  by  a  recently  published  review  of 
drug  metabolism  in  lungs.  (Brandenburger  Brown,  E.A.,  The  localization,  metabolism, 
and  effects  of  drugs  and  toxicants  In  lung.  Drug  Metabolism  Reviews,  3,  33,  1974; 
Chiou,  W.L.,  Potential  pitfalls  in  the  conventional  pharmaco- 
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kinetic  studies:  effects  of  the  initial  mixing  of  drug  In  blood  and  the  pulmonary 
first-pass  elimination.  Journal  of  Pharmacokinetics  and  Biopharmaceutics,  7,  527, 
1979;  Welch,  R.M.,  Cavallito,  J.,  Loh,  A.,  Effect  of  exposure  to  cigarette  smoke 
on  the  metabolism  of  benzo (a) pyrene  and  acetophenetidin  by  lung  and  intestine  of 
rats,  Toxicol.  Appl.  Pharmacol.,  23,  749,  1972).  The  second  possibility  is  also 
supported  by  the  literature  (Fiserova-Bergerova ,  V.,  Teisinger,  J.,  Pulmonary 
styrene  vapor  retention.  Industrial  Medicine  and  Surgery,  34,  620,  1965;  Bardodej , 
Z.,  The  value  and  use  of  exposure  tests  XIL.  Mercury  exposure  test,  Ceskoslovenska 
Hygiena,  8,  157,  1963;  Landahl,  H.D.,  Hermann,  R.G.,  Retention  of  vapors  and  gases 
in  the  human  nose  and  lungs.  Archives  of  Industrial  Hygiene  and  Occupational  Medi¬ 
cine,  1,  36,  1950).  However,  quantitative  description  of  adsorbtlon  and  desorbtion 
in  respiratory  airways  is  lacking.  We  were  attempting  to  fill  this  gap.  The  third 
possibility,  the  limitation  of  uptake  rate  by  diffusion,  is  supported  by  older 
literature  (Kety,  S.S.,  The  theory  and  applications  of  the  exchange  of  inert  gas  at 
the  lungs  and  tissues.  Pharmacological  Reviews,  3,  1,  1951;  Hunter,  A.R.,  The  group 
pharmacology  of  anaesthetic  agents.  I:The  adsorbtlon-el imination  of  Inhaled 
drugs,  Brit. J.Anaesth. ,  28,  244,  1956).  However,  diffusion  Is  not  our  priority 
concern,  since  there  was  no  systematic  deviation  of  the  ratio  of  concentrations 
In  arterial  blood  and  alveolar  air  from  blood-air  partition  coefficient. 

3)  Interference  of  simultaneously  administered  xenobiotics.  We  ^^rved  that 
if  two  vapors  are  Inhaled  simultaneously,  the  one  reduces  pulmonary  uptake  of  the 
other.  We  expressed  the  hypothesis  that  reduced  uptake  Is  caused  by  competitive 
Inhibition  of  metabolism.  Experiments  corroborating  this  hypothesis  are  underway, 
but  so  far  no  data  suitable  to  application  for  modeling  was  obtained. 

Significance: 

With  progressing  technology,  toxicologists  and  hygienists  confront  the  problem 
of  securing  safe  exposure  to  an  increasing  number  of  air  pollutants.  The  adverse 
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biological  effect  of  pollutants,  like  the  therapeutic  effect  of  drugs,  is  related 
to  blood  concentration  and/or  time  integral  of  concentration  in  the  target  organs. 
Passage  of  pollutants  from  the  environment  to  the  target  organ  has  the  same  signifi¬ 
cance  as  migration  of  drugs  from  site  of  administration  to  the  target  organ.  Pharma¬ 
cokinetics  describing  the  transport  of  drugs  in  the  body  is  a  potent  tool  for  design¬ 
ing  dosage  regimens  of  optimum  therapeutic  effect.  We  employed  similar  methods  to 
propose  exposure  regimens  with  equivalent  adverse  effect.  The  evaluation  is  based 
on  comparison  of  concentrations  reached  in  arterial  blood.  Two  of  our  observations 
deserve  special  attention:  1)  The  first  observation  was  made  while  using  the  program 
for  simulation  of  short-term  increase  of  exposure  concentration.  It  was  calculated 
that  the  rising  of  vapor  concentration  in  tissues  depends  on  the  solubility  and  eli¬ 
mination  of  the  substance.  It  was  demonstrated  that  if,  during  the  first  four  hours 
of  an  eight  hour  exposure,  the  exposure  concentration  rises  for  a  short  period  (15 
minutes)  two  or  three  times,  tissue  concentrations  of  well  soluble  vapors  at  the 
end  of  excursion  are  smaller  than  at  the  end  of  exposure.  Based  on  this  simulation, 
we  suggested  to  the  Threshold  Limited  Value  Airborne  Contaminants  Committee  of 
American  Conference  for  Industrial  Hygienists  that  short-term  excursion  limits  for 
industrial  pollutants  (STEL)  should  be  based  on  solubility  of  the  vapor,  rather 
than  on  the  average  concentration  (TWA),  as  is  the  current  practice.  Similar  sugges¬ 
tion  is  made  by  Japanese  investiaators  (Koizumi, A.,  Sekiguchi,  T.,  Konno,  M.,  Ikeda, 

M. ,  Evaluation  of  the  time  weighted  averaqe  of  air  contaminants  with  sDecial  referen¬ 
ces  to  concentration  fluctuation  and  biological  half-time,  American  Industrial  Hygiene 
Association  Journal,  41,  693,  1980). 

2)  The  second  important  observation-  the  nonlinear  dependence  of  biological 
effect  on  exposure  concentration -concerns  the  extrapolation  of  toxicological  data. 

The  observed  limited  capacity  of  elimination  process  can  affect  the  elimination  not 
only  quantitatively  but  also  qualitatively.  Gehrinq  and  Blau  discussed  mechanism 
of  carcinogenesis  with  regard  to  dose  dependent  changes  in  metabolic  pathways 
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of  inhaled  carcinogens  (Gehring,  P.J.,  Blau,  G.E.,  Mechanisms  of  carcinogenesis: 
dose  response.  Journal  of  Environmental  Pathology  and  Toxicology.  1,  163,  1977). 

We  are  showing  In  our  experiments  with  halothane  that  similar  changes  in  metabolic 
pathways  can  be  induced  if  the  vapor  is  inhaled  simultaneously  with  another  vapor. 
The  studies  show  that,  in  addition  to  enhancing  secondary  elimination  pathways,  the 
saturation  of  the  major  elimination  pathways  results  In  the  rapid  rising  of  vapor 
concentration  in  tissues.  The  quantitative  and  qualitative  changes  introduced  by 
capacity  limited  processes  render  dubious  the  extrapolation  of  safe  exposure  concen¬ 
tration  from  toxicological  studies  performed  at  higher  exposure  concentrations  or 
from  studies  in  which  other  than  inhalation  administration  was  used.  In  order  to 
prevent  an  unexpected  reaction  in  the  presence  of  additional  xenobiotics  (Includ¬ 
ing  drugs)  exposure  concentrations  below  those  which  saturate  detoxifying  mechanisms 
should  be  recommended  whenever  determining  the  safety  limit  for  occupational  or 
environmental  exposure. 

It  is  highly  desirable  to  continue  to  study  the  parameters  defining  nonlinear 
elimination,  and  to  prepare  a  eneral  simulation  model  which  accomodates  capacity 
limited  processes. 


* 


Program  LINEAR  (also  POMOC) 
General  Linear  model  for  5  compartments 


N,NGX 
gci>,c(I),i=i...n 
6X<n,I»l...NGX 
(  IF  NGX=0 ,  NO  GX(I)  IS  RERD) 
STEP 

NUMBER  OF  TIME  SECTIONS, NUMSEC 
E(I),TIMENO<  I  >,!=!..  .NUMSEC 


52  ORTfl  44,128 

53  ORTB  28.3,688 

54  DRTR  11.2,5000  TABLE  1 

58  DRTR  6,12 

60  DRTR  8.47 
62  DRTR  480 
64  DRTR  2 
66  DRTR  1,488 
68  DRTR  0,960 

100  DIM  R(  6,& ),B< 6,6,7 ),C( 6  ),CC( 7 ) 

110  DIM  G<6),GX(6),XT<6),XVNTC6),CURC6) 
120  DIM  CVNTC3),TD(30),MSTEC30),E<30> 
130  DIM  X(6),RR(6),XM<6) 

135  DIH  W<6) 

140  DIM  BB(7),H(7),D<  11), 2(200) 

150  DIM  TSC 100 ),NRP( 100,6) 

160  OIM  PROUD( 100,6 ),MZ( 100,8) 

170-  RERD  N,NGX 

180  PRINT  “N=“,N 

181  PRINT  "NGX=",NGX 
190  N1  *  N  +  1 

200  N2  *  N  -  1 
210  N5  *  NGX  +  1 
220  NS  =  NGX  +  2 

230  N7  *  NGX  +  3 

231  NM1  »  N  -  1 


240 

■aPMKitiEVfc 

250 

FOR  I  a  1  TO  N 

260 

XM<  I  )  *  0 

270 

XVNT(I)  a  0 

280 

X(I)  a  0 

285 

M(I)  =  0 

290 

GX(  I  )  =  0 

300 

NEXT  I 

305 

PRINT 

LINEAR 

2 


310 

rcm  pr  on  ci  cmcmtr 

311 

PRINT  "G< I > 

320 

FOR  I  =  1  TO  N 

330 

REPO  Gfl),ca> 

348 

PRINT  G< I >,C< I ) 

350 

NEXT 

360 

IF  NGX  *  0  GOTO  420 

361 

PRINT 

362 

PRINT  "NGXC I )" 

370 

FOR  I  a  1  TO  NGX 

380 

READ  GX<I) 

390 

PRINT  I,GX<n 

400 

NEXT  I 

402 

RFM  START  PRINTER . 

404  INPUT  "WANT  PRINTER?  <Y/NVjR* 

406  IF  R*  *  "V"  THEN  PR#  1 
416  RFM  CaLQlgEE  OSVHPTQTS 
420  GOSUB  5000 

425  W<N>  *  G<N> 

426  PRINT 

430  GOSUB  5500 

431  PRINT  "ASYMPTOTIC  UPLUES  FOR  UNIT 

432  PRINT  H .  CONCENTRATION" 

440  FOR  I  *  1  TO  N 

450  PRINT  I,X<I> 

460  NEXT  I 
470  FOR  I  =  1  TO  N 
480  X<I>  *  0 
490  W<I>  *  0 
510  NEXT  I 
610  PRINT 

615  REH  SIMULATION  ..IMFQRMflTIQfcL 
620  READ  SEP 
630  PRINT  "STEP-",SEP 
640  READ  NUS 

650  PRINT  "#  OF  TIME  SECTIONS**", NUS 

654  PRINT 

655  PRINT  "CONCENTR.  TIME  IN  MIN" 

660  FOR  I  »  1  TO  NUS 

670  repo  Ecn/rocn 

686  PRINT  E( I >,TD< I > 

690  NEXT  I 

700  FAN  *  G(N)  /  C<N5 

710  GOSUB  3000 

720  GOSUB  3160 

730  GOSUB  3890 

740  T  a  0 

750  IT  a  0 

760  UOLO  a  0 

770  F2  a  0 

780  FOR  I  *  t  TO  NUS 
790  MSTE<  I)  =  TD<n  /  SEP 
800  NEXT  I 
810  PRINT 


LINEAR 

3 


<•  * 


320  REM  STPPT  nr  PTHtjt  QTTON 
573  r^R  JJ  =  1  TO  si  ts 

340  MM  =  MSTECJJ) 

350  FOT  =  EC JJ) 

360  FOR  JK  a  1  TO  MM 
370  F7  a  F2  +  FOT  *  SEP 
380  T  s*  T  +  SEP 
890  IT  »  IT  +  1 
900  TSCIT)  =  T 

910  REM  CALCULATE  XCI),XVNTCT)  -  CQNCENTRATT  0M5 

920  FOR  L  =  1  TON 

330  SI  =8 

940  S2  *  0 

950  FOR  M  =  1  TO  N 

960  SI  *  SI  +  RCL#M)  *  XCM) 

970  S2  *  S2  +  BCL,M,N1)  *  XCM) 

980  NEXT  M 

990  XTCL)  =  SI  +  CCL)  *  FOT 
1000  U  =  S2  +  CCCL)  *  FOT 
1010  XVNTCL )  *  IJ  +  XMCL) 

1020  NEXT  L 

1030  REM  TRANSFER  XC I ),XYNTC I ) . 

1040  FOR  L  =  1  TO  N 
1050  XMCL)  a  XVNTCL) 

1060  XCL)  =  XTCL) 

1070  NPP<IT,L)  =  XCL) 

1060  NEXT  L 

1030  REM  CALDJLATE  CURC  I ) ,CYNTC  T  )  -  ROTES  AND  AMOUNTS 
1100 ic  N6X  =  0  THEN  $0*0  11?U 
1110  FOR  r  =  1  TO  NGX 
1120  CURC  I)  =  GXC  I  )  *  XCI) 

1130  PROUDC IT# I )  =  CURC I) 

1140  CYNTCI)  a  GXC  I)  *  XVNTCI) 

1150  WZCIT,I)  =  CVNTCI) 

1160  NEXT  I 

1170  CURCN5)  a  GCN)  *  CFOT  -  XCN)) 

1180  PROUDC IT,N5)  a  CURCN5) 

1190  CVNTCN5)  *  GCN)  *  CF2  -  XVNTCN)) 

1200  N2CIT,N5)  =  CVNTCN5) 

1210  U  =  CYNTCN5)  -  UOLD 
1220  IF  U  >  0  THEN  CVNTCN6)  =  CVNTCN6)  +  U 
1230  IF  U  <  0  THEN  CVNTCN7)  =  CYNTCN7)  +  U 
1240  UOLD  =  CVNTCN5) 

1250  WZC IT *N6 )  a  CVNTCN6) 

1260  WZC  IT*N7)  a  CVNTCN7) 

1265  REM  DISPLAY  DATA  CSEE  OUTPUT  TABLE) 


TABLE  2 


1269  PRINT 

1275"' 'WYRt  HTIM^a»T 

1280  PRINT  XC  1  )#XC2) 

1281  PRINT  XC3),XC4) 
1290  PRINT  CURC 1 )»CUR<2) 
,1300  PRINT  CVNTCI) 


1320  NEXT  JK 
1330  NEXT  JJ 


C  4) 


LINEAR 

4 


1 


I 


1225...  Bgt  FILS  EftTfi 

1339  PRINT 

1340  INPUT  "WONT  FILE*?1  <V^N)";Rf 


1345  IF  Rt  a  "N"  THEN  END 
1350  Ff  a  “RESULT" 

1360  0*  «  CWRt  (4) 

1370  PRINT  Df;"OPEN";Ft 
1380  PRINT  Of; "DELETE" ;Ff 
1390  PRINT  Of;"OPEN";Ff 
1400  PRINT  Of ; "WRITE" ;Ff 
1410  PRINT  N 
1420  PRINT  N5 
1430  PRINT  N7 
1440  PRINT  IT 
1450  FOR  I  a  1  TO  IT 
1460  PRINT  TS< I > 

1470  FOR  J  *  1  TO  N 
1480  PRINT  NflP(I,J) 

1490  NEXT  J 

1500  FOR  J  *  l  TO  N5 

1510  PRINT  PROUD(I,J) 

1520  NEXT  J 
1530  FOR  J  *  1  TO  N7 
1540  PRINT  W2< I 
1550  NEXT  J 
1560  NEXT  I 

1570  PRINT  Of; "CLOSE" ;Ff 
1575.  REH  ORGQNISE  DflTQ 
1590  Z<0)  a  IT 
1600  2(1 )  =  0 
1610  FOR  JM  a  1  TO  3 
1620  IF  JM  =  1  THEN  NK  =  N 
1630  IF  JM  =  2  THEN  NK  a  N5 
1640  IF  JM  =  3  THEN  NK  =  N7 
1650  FOR  J  *  l  TO  NK 
I860  IF  JM  a  1  THEN  GOSUB  4670 
1670  IF  JM  a  2  THEN  GOSUB  4740 
1680  IF  JM  =  3  THEN  GOSUB  4810 
1690  FOR  I  =  1  TO  IT 
1700  MM  a  2  *  I 
1710  Z(MM)  =  TS(I) 

1720  MM  a  MM  +  1 

1730  IF  JM  a  l  THEN  ZCMM)  =  NflP(I,J) 

1740  IF  JM  a  2  THEN  2<MM)  a  PROUDC I#J5 

1750  IF  JM  a  3  THEN  2(MM)  =  WZ(I.J) 

1760  NEXT  I 

1770  GOSUB  4560 

1780  NEXT  J 

1790  NEXT  JM 

1800  END 


L 


RtM  rfjOMtJI  QTTpij  pc  Tt-IC  CypTO-j  MQTp 


3560  LIH  *  100 
3570  GOSUB  4120 


3589  REH  CALCULATE  THE  RESIDUES 

3590  FOR  L  =  1  TO  N 
3600  RL  =  RRCL) 

3610  PROD  =  1 

3620  FOR  K  =  1  TO  N 

3630  IF  L  =  K  GOTO  3650 

3640  PROD  =  PROD  *  (PL  -  RRC K  )  ) 

3650  NEXT  K 

3660  CCL)  =  PROD 

3670  NEXT  L 

3680  FOR  I  =  1  TO  N 

3690  FOR  J  =  1  TO  N 

3700  FOR  K  =  1  TO  N 

3710  COCK)  =  BCI,J,K) 

3720  NEXT  K 
3730  FOR  L  *  1  TO  N 
3740  RL  «  RRCL) 

3750  PUB  ■  CCC1) 

3760  FOR  K  =  2  TO  N 
3770  PUR  a  PUB  *  RL  +  CCCK) 

3780  NEXT  K 


3810 

3820 

3825 

3826 
3830 


TIME  CONST." 


3790  BCI,J,L)  a  PUR  /  CCL) 

3800  NEXT  L 
NEXT  J 
NEXT  I 
PRINT 

PRINT  "EXPONENT 
FOR  L  =  1  TO  N 
3835  U  =  -  1  ✓  RRCL) 

3840  PRINT  RRC  L )*U 
3850  FOR  I  =  1  TO  N 
3860  NEXT  I 
3870  NEXT  L 
3880  RETURN 

~180  REN  CPLCULRTE  EXPONENT I RLS  RND  IN' 
3900  FOR  I  =  1  TO  N 
3910  CCCI)  a  0 


3920  CCI)  a  0 
3930  FOR  J  =  1  TO  N 
3940  PCI,J)  =  0 
3950  BCI,J,N1)  =  0 
3960  NEXT  J 
3970  NEXT  I 


3980  FOR  K  =  1  TO  N 

3990  EU  a  EXP  <RR<fO  *  SEP) 

4000  EU  =  CEU  -  1)  /  RRCK) 

4010  EH  »  (EU  -  SEP)  ✓  RRCK) 

4020  FOR  I  =  1  TO  N 

4030  CCI)  a  CCI)  *  BCI,N,K)  *  EU  »  FPN 

4040  can  =  OOC  I  t  +  P(I,N,K)  *  Ell  *  FPN 

4050  FOR  J  a  i  JO  H 

4Rp0  ccr,j’»  =  +  bci,j,k)  *  £■  ■ 

4070  BCI,J,N1)  a  BCI,J,N1)  +  3CI^,K)  *  Pit 

4080  NEXT  J 

4090  NEXT  I 

4100  NEXT  K 

4110  RETURN 


_ REM  ROOT-FTNPIMg  STEERING  routine 

4130  for  i  *  i  Toni 

4140  BBC  I  )  =  CCCI) 

4150  NEXT  I 
4160  XX  =  0 
4170  NO  =  N 
41S0  NP  =  Ml 

4190  FOR  IROOT  =  1  TO  N 

4200  ITER  =  0 

4210  GOSUB  4420 

4220  ITER  =  ITER  +  1 

4230  IF  ITER  c  LIM  SOTO  4260 

4240  PRINT  "ITERATIONS  EXCEEDED" 

4250  END 

4260  OX  »  -  FU  /  ON 
4270  XX  =  XX  +  DX 
4280  U  =  BBS  COX) 

4290  IF  U  >  EPS  THEN  SOTO  4210 

4300  RRC IROOT)  =  XX 

4310  GOSUB  4360 

4320  NO  *  NO  -  1 

4330  NP  a  NP  -  l 

4340  NEXT  IRROT 

4350  RETURN 


4370  IF  NP  <  2  THEN  GOTO  4410 

4380  FOR  K  =  2  TO  NP 

4390  BBCK)  =  BBCK)  +  BBCK  -  1 )  *  XX 

4400  NEXT  K 

4410  RETURN 

4420  REM  POLVNOMIRL  EURLUATION 
4430  FOR  I  =  1  TO  NP 
4440  HC I )  a  BBC  I ) 

4450  NEXT  I 

4460  FOR  I  =  2  TO  NP 

4470  HCI)  a  HCI)  +  XX  *  H<  I  -  1) 

4480  NEXT  I 

4490  FU  a  HCNP) 

4500  IF  NO  =  1  THEN  GOTO  4540 

4510  FOR  I  =  2  TO  NO 

4520  HCI)  a  HCI)  +  XX  *  HCI  -  1  ) 

4530  NEXT  I 

4540  ON  a  HCNO) 

4550  RETURN 


3  f-UPt  ;•  i 

4570 

’  PRINT  Dfi "OPEN" iFf 

4580 

PRINT  Dfi "DELETE" iFf 

4590 

PRINT  Dfi "OPEN" iFf 

4800 

PRINT  Dfi "WRITE" ;Ff 

4610 

ML  =  2  *  IT  +  1 

4620 

FOR  I  =  0  TO  ML 

4630 

PRINT  2( I ) 

4649 

NEXT  I 

4658 

PRINT  Dfi "CLOSE" iFf 

4660 

RETURN 

4670 

IF  J  =  1  THEN  Ff  = 

"NflPl" 

4680 

IF  J  =  2  THEN  Ft  = 

"NRP2" 

4690 

IF  J  =  3  THEN  Ft  = 

"NRP3" 

4700 

IF  J  =  4  THEN  Ff  = 

"NRP4" 

4710 

IF  J  =  5  THEN  Ff  a 

"NRP5* 

4720 

IF  J  =  6  THEN  Ft  = 

"NRP6" 

4730 

RETURN 

4740 

IF  J  a  1  THEN  Ff  = 

"PR0UD1" 

4750 

IF  J  a  2  THEN  Ft  = 

"PR0UD2" 

4760 

IF  J  =  3  THEN  Ff  = 

"PR0U03* 

4770 

IF  J  =  4  THEN  Ff  = 

"PR0U04" 

4780 

IF  J  a  5  THEN  Ft  a 

"PR0UD5" 

4790 

IF  J  a  6  THEN  Ft  = 

"PR0UD6" 

4800 

RETURN 

4810 

IF  J  a  1  THEN  Ff  a 

"INGR1" 

4820 

IF  J  a  2  THEN  Ff  = 

"INGR2" 

4830 

IF  J  a  3  THEN  Ff  * 

"INGR3 

4840 

IF  J  a  4  THEN  Ft  a 

"INGR4" 

4850 

IF  J  =  5  THEN  Ff  a 

"INGR5" 

4860 

IF  J  a  6  THEN  Ft  = 

" INGR6" 

4870 

IF  J  a  7  THEN  Ft  a 

"INGR7" 

4880 

IF  J  a  8  THEN  Ft  = 

"INGR8" 

4890 

RETURN 

5000  REM  LU  FQCT0RI2QTTQN 
5010  D(N)  =  0 
5020  FOR  I  =  1  TO  NM1 
5020  D<I)  =  G(I)  +  GXCI) 

5040  L  »  N  +  I 
5050  D(L)  =  -  G<  I  >  /  DC  I) 

5060  O(N)  ■  O(N)  +  G(I)  *  (1  +  O(L)) 

5070  NEXT  I 

5080  OCH)  =  0(N )  +  6(N) 

5090  RETURN 

,5500  REM  SUBROUTINE  LUSOL 
5520  2(N)  =  W(N> 

5530  FOR  I  »  1  TO  NH1 

5535  2(N)  =  2(H)  *  2(1)  *  G<  I) 

5540  2(1)  *  W(I)  /  0(1) 

5550  NEXT  I 

5560  2(N)  *  2(N)  /  0(N) 

5570  U  =  2(N) 

5580  FOR  I  »  1  TO  NM1 

5590  L  a  n  +  I 

5600  X(I)  *  2(1)  -  0(L)  *  U 

5810  NEXT  I 

5629  X(N)  s  !_! 
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DETERMINATION  OF  KINETIC  CONSTANTS 
FROM  PULMONARY  UPTAKE 

V.  Fiserova-Bergerova,  Ph.D. 

University  of  Miami  School  of  Medicine 
Miami ,  Florida 


The  adverse  biological  effect  of  an  air  pollutant  is  related  to  its 
concentration  and/or  time  integral  of  concentration  in  the  target  organ. 
The  passage  of  pollutant  from  the  environment  to  the  target  organ  is 
a  dynamic  process,  determined  by  the  physical  and  chemical  properties 
of  the  pollutant  and  by  the  physiological  parameters  of  the  exposed 
subject  (Eger,  1963).  Inhaled,  nonwater  soluble  vapors  are  removed 
from  the  body  by  pulmonary  and  metabolic  clearance.  Since  pulmonary 
clearance  takes  place  only  after  the  end  of  exposure  (or  when  exposure 
concentration  decreases),  the  vapor  is  removed  during  exposure  only 
by  metabolic  clearance  (Fiserova-Bergerova  et  ai. ,  1974). 

Metabolic  clearance  in  vivo  is  usually  defined  by  the  half-time  of 
disappearance  of  the  xenobiotic  from  plasma  after  bolus  administration, 
or  as  the  half-time  of  urinary  excretion  of  xenobiotic  metabolites.  The 
determination  of  plasma  clearance  requires  frequent  blood  sampling, 
which  imposes  stress  on  subjects.  The  urinary  excretion  of  metabolites 
is  affected  by  a  variety  of  factors,  such  as  distribution  of  metabolites 
in  the  body,  binding,  and  renal  clearance. 

Inhalation  administration  makes  possible  a  noninvasive,  accurate 
measurement  of  the  rate  of  overall  metabolism  of  inhaled  vapors 
(Teisinger  and  Soucek,  1952). 

The  pulmonary  uptake  rate,  u,  is  the  sum  of  retention  rates  of 
vapor  in  tissues,  utjs,  and  the  rate  of  overall  metabolism,  um: 


+  u 


m 


(1) 


! 

I 


The  pulmonary  uptake  can  be  determined  from  the  difference 
between  vapor  concentrations  in  inhaled  and  mixed- exhaled  air  multiplied 
by  minute  ventilation: 


u  =  (C 


inh 


Cexh) 


(2) 
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where  Cjnh  and  CeXh  are  vapor  concentrations  (mg /liter)  in  inhaled  air 
and  in  mixed -exhaled  air.  and  V  is  minute  ventilation  (liter /minute) . 


We  have  presented  a  compartments]  model  which  can  be  used  to 
determine  the  retention  rate  of  inhaled  vapor  in  tissues.  Using  three 
comp.artments  (VRG  -  vessel  rich  tissues.  MG  -  muscles  and  skin,  FG  - 
fat  and  fat  marrow),  the  retention  rate: 


u 


tis 


FVRGXbl/air  exp 


(• 


VRG 

T 


VRG  VRG /bl 


+  FMGXbl/air  SXp 
+  FFGXbl/air  exp 


(■ 

(- 


fmc  \ 

V  T  4  I 

VMG  MG/bl  / 

FFG  \ 

V  X  1  j 

VFG  FG/bl  /  j 


(3) 


where  F's  are  blood  flows  (liter /minute)  through  the  compartments,  V's 
(1)  are  their  volume.  .X's  are  the  corresponding  partition  coefficients  of 
inhaled  vapor  at  37°C,  Cajv  is  the  vapor  concentration  (mg /liter)  in 
alveolar  air,  and  exp  is  the  natural  logarithm. 

Substituting  from  Equations  2  and  3  in  Equation  1,  the  rate  of' 

overall  metabolism  can  be  determined.  Since  the  retention  of  vapor  in 

tissues  is  an  exponential  function,  the  retention  rate  diminishes  with 

exposure  duration.  This  determination  is  most  accurate  during  apparent 

steady  state,  when  retention  by  tissues  is  small  (ii  >>  u  .  ) . 

m  tis 

We  demonstrated  the  effect  of  metabolism  on  pulmonary  uptake  in 
an  informed  volunteer  patient  who  was  anesthetized  with  fluroxene 
(CH2-CH.O.CH 2.CF3)  in  the  presence  of  a  small  concentration  of  non- 
metabolized  isoflurane  (CHFZ. )  .CHC1.CF  3)  (Fiserova-Bergerova  and 
Holaday ,  1979).  The  amount  of  metabolites  accounts  for  less  than  1% 
of  isoflurane  uptake  (Holaday  et  al.,  1975),  and  about  45%  of  fluroxene 
uptake  (Gion  et  al.  ,  1974). 

During  anesthesia,  samples  of  inhaled  gas  and  end-exhaled  gas 
were  drawn  during  the  appropriate  phase  of  respiration  via  a  nylon 
cannula  inserted  in  the  endotracheal  tube  at  4  to  10  minute  intervals. 
Mixed-exhaled  gas  was  obtained  at  the  same  time  at  the  outlet  of  a 
mixing  chamber  interposed  in  the  expiratory  breathing  tube.  Gas 
samples  were  collected  in  20  ml  glass  syringes.  At  the  same  time. 
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f  luroxene 


minute  ventilation  was  measured  with  a  Wright  respirometer.  The  cumu¬ 
lative  uptake  0  was  determined:  (1)  from  the  amount  of  anesthetics 
delivered  by  syringe  in  the  closed  anesthetic  circuit,  and  (2)  from  the 
sum  of  the  differences  between  concentrations  of  inhaled  and  mixed- 
exhaled  air  multiplied  by  minute  ventilation,  and  time  intervals  between 
sampling  (t  in  minutes) 


D 


E  V  t  (C. 


inh 


-  C  .  ) 
exh 


(4) 


The  cumulative  uptake  predicted  by  integration  of  Equation  3  as  reten¬ 
tion  in  tissues  correlates  with  the  measured  uptake  of  isoflurane  (calcu¬ 
lated  by  Equation  4) ,  but  the  measured  fluroxene  uptake  greatly  exceeds 
the  calculated  fluroxene  retention  in  tissues  (Figure  1).  The  difference 
accounts  for  fluroxene  metabolism. 


Figure  1.  Cumulative  uptake  of  isoflurane  (i)  and  fluroxene  (f)  administered 
simultaneously  to  a  surgical  patient.  Cumulative  uptake  is  plotted  against  the 
time  after  the  start  of  anesthesia.  The  dashed  lines  are  uptake  curves  pre¬ 
dicted  by  integration  of  Equation  3  for  alveolar  concentrations  of  fluroxene 
(154  mg/liter)  and  isoflurane  (3  mg/liter). 


To  study  the  effect  of  exposure  concentration  on  pulmonary 
uptake,  the  following  assumptions  were  made:  (1)  the  retention  of 
vapor  in  tissue  is  a  first  order  process,  which  means  its  rate  constant 
is  concentration  independent:  and  (2)  metabolic  clearance  is  a 
limited-capacity  process  described  by  Michaelis-Menten  kinetics. 

To  determine  the  Michaelis-Menten  constants  in  vivo,  we  exposed 
male  rhesus  monkeys  (approximately  3  kg)  consecutively  to  three  con¬ 
centrations  of  one  of  the  following  compounds:  benzene,  halothane 
(CF  j-CHClBr) ,  methylene  chloride,  or  trichloroethylene.  Concentrations 
were  in  the  range  of  TLV  during  the  first  exposure;  equal  to  five  times 
TLV  during  the  second  exposure;  and  equal  to  25  TLV  during  the;  third 
exposure.  In  order  to  reach  apparent  steady  state,  each  exposure^  lasted 
approximately  two  and  one- half  hours.  Vapors  were  administered  in  light 
sernylane  anesthesia  via  endotracheal  tube.  The  following  parameters 
were  measured:  (1)  vapor  concentrations  in  inhaled  air  (Cexp)»  mixed- 
exhaled  air  (Cexh)>  and  end-exhaled  (Calv)  and  arterial  blood(Cart) ; 

(2)  blood-gas  partition  coefficients  (Xbl/air);  (3)  minute  ventilation; 

(4)  blood  pressure  and  pulse  rate;  and  (5)  blood  gases  and  PC02  in 
mixed-exhaled  air.  Uptake  rate,  metabolic  rate,  and  alveolar  ventilation 
were  calculated  from  the  measured  data.  Apparent  Michaelis-Menten 
constants  of  overall  metabolism  in  vivo  (Km)  were  calculated  from  double 
reciprocal  plots  of  metabolic  rate  versus  Calv,  Cexp»  or  C^xt/ Xbl/air 
(measured  at  a  steady  state) ,  and  versus  calculated  concentrations  in 
tissues. 

In  Figure  2,  the  double  reciprocal  plots  from  benzene  and  methylene 
chloride  are  presented.  Km  values  related  to  the  concentrations  in  alveo¬ 
lar  air  at  steady  state  for  all  four  studied  compounds  are  in  Table  1. 

Figures  3  and  4  demonstrate  double  reciprocal  plots  of  metabolic 
rates  of  trichloroethylene  and  halothane  versus  exposure  concentra¬ 
tion,  alveolar  concentration,  arterial  concentration,  and  tissue  concen¬ 
tration.  The  data  indicate  that  Km  values  depend  on  the  site  in  which 
the  concentration  is  measured,  but  the  maximum  metabolic  rate  (Vmax) 
is  the  same  regardless  of  whether  it  is  derived  from  concentration  in 
inhaled  air,  alveolar  air,  arterial  blood,  or  tissue. 


Figure  2.  Double  reciprocal 
plot  of  uptake  rates  of 
benzene  and  methylene 
chloride  versus  alveolar 
concentration.  The  lines 
represent  optimum  fit  to 
experimental  data  obtained 
in  rhesus  monkeys  (3  kg 
males) . 


TABLE  1.  APPARENT  MICHAEUS-MENTEN  CONSTANTS  IN  VIVO 
(3  KG  MALE  RHESUS  MONKEY) 


Compound 

Vmax  (mg/rain) 

Km*  (mg/liter) 

TLV**  (mg /liter) 

Benzene 

0.3 

0.6 

0.03 

Trichloroethylene 

0.7 

1.0 

0.27 

M  Jthylene  Chloride 

0.7 

1.1 

0.36 

Halothane 

1.2 

7.4 

0.40 

♦Related  to  Caiv 

•♦Threshold  limit  values  recommended  by  ACGIH  in  1979  as  TWA. 
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Figure  3.  Double  reciprocal  plot  of  uptake  rates  of  trichloroethylene 
versus  measured  concentrations  in  inhaled  air  (CeXp),  end-exhaled 
air  (Calv)*  arterial  blood  (Cart) .  and  calculated  concentrations  in 
tissues . 
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The  lines  represent  optimum  fit  to  experimental  data  obtained  in  rhesus 
monkeys  (3  kg  males). 
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Figure  4.  Double  reciprocal  plot  of  uptake  rates  of  halothane  versus 
measured  concentrations  in  inhaled  air  (Cexp) ,  end-exhaled  air  (Calv) 
arterial  blood  (Cart)*  and  calculated  concentrations  in  tissues. 
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The  lines  represent  optimum  fit  to  experimental  data  obtained  in  rhesus 
monkeys  (3  kg  males). 


The  differences  in  Km  values  might  be  explained  by  the  three 
concentration  gradients  on  the  pathway  of  the  vapor  from  the  en¬ 
vironment  to  the  metabolic  site:  (1)  Vapor  concentration  entering 
the  lung  with  each  breath  is  smaller  than  the  exposure  concentration 
because  of  the  dilution  by  alveolar  air  from  deadspace;  (2)  When 
the  air  reaches  the  alveoli,  the  partial  pressures  in  alveolar  air  and 
arterial  blood  are  readily  equilibrated.  Uptake  of  vapor  by  arterial 
blood  reduces  the  concentration  further.  The  concentration  decreases, 
depending  on  the  cardiac  output  and  alveolar  ventilation,  on  the  solu¬ 
bility  of  vapor  in  blood,  and  on  the  concentration  of  vapor  in  mixed- 
venous  blood;  (3)  Arterial  blood  transfers  the  vapor  to  the  tissues, 
where  it  is  retained  to  the  extent  that  partial  pressures  in  tissue 
and  venous  blood  are  equilibrated.  The  concentration  gradient  • 
Cart/Cven  *n  ^lood  which  supplies  metabolic  sites  is  further  increased 
by  metabolic  clearance. 

At  steady  state,  the  partial  pressures  of  nonmetabolized  vapor 
equilibrate,  and  the  vapor  concentrations  in  tissues  equal  the  exposure 
concentration  multiplied  by  the  appropriate  partition  coefficient.  If  the 
vapor  is  metabolized,  the  concentrations  are  reduced. 

Employing  our  nonlinear  model  (Fiserova-Bergerqva  et  al. ,  in 
preparation) ,  we  examined  the  conditions  under  which  this  method  for 
Kiq  and  Vmax  is  applicable,  and  found  the  following  limitations: 

(1)  Metabolism  must  be  concentration  dependent.  According  to 
Michaelis-Menten  kinetics,  this  requirement  is  met  if  substrate  concen¬ 
trations  are  smaller  than  10  Km-  This  means  that  the  studies  must  be 
performed  in  the  range  of  exposure  concentrations  which  are  smaller 
than  10  Km . 


C 


exp 


<  10  K 
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(2)  The  system  cannot  be  flow-limited.  •  This  requires  that  trans¬ 
portation  rate  of  vapor  from  environment  to  the  metabolic  site  is  larger 
than  metabolic  rate.  This  condition  is  met  if: 


V 

max 


F  Abl/air  if  , 
_  alv 

+  F  Abl/air 

Valv 


This  expression  can  be  rearranged: 
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Optimum  conditions  for  determination  are:  (a)  the  vapor  is  highly 
susceptible  to  biotransformation,  and  (b)  the  metabolite  sites  are  well 
perfused  or  the  vapor  is  well  soluble  in  blood. 

We  analyzed  the  tissues  of  rats  and  monkeys  exposed  to  different 
concentrations  of  halothane  and  trichloroethylene  to  determine  the  effect 
of  exposure  concentrations  on  concentrations  of  these  vapors  in  tissues. 
When  exposure  concentrations  were  larger  than  Kn,  the  metabolic 
clearance  diminished  and  the  concentration  ratios  Cexp/Ctis  increased 
(Fiserova-Bergerova,  unpublished  data).  The  same  conclusions  were 
drawn  using  our  nonlinear  mathematical  model  (Fiserova-Bergerova 
et  al.  .  in  preparation) . 


The  determination  of  metabolic  rate  from  pulmonary  uptake  is  also 
suitable  for  studying  the  effect  of  modifiers  on  metabolism  of  inhaled 
vapors.  In  experiments  similar  to  those  described  above,  we  adminis¬ 
tered  two  vapors  simultaneously  to’  monkeys.  One  vapor  was  adminis¬ 
tered  at  a  low  constant  concentration:  the  concentration  of  the  other 
vapor  -  the  modifier  -  was  increased  in  three  steps.  Data  from  these 
experiments  are  in  Figure  5.  As  the  concentration  of  modifier  increased, 
the  pulmonary  uptake  of  the  studied  vapor  decreased.  The  decrease  is 
probably  caused  by  competitive  inhibition  of  metabolism  of  the  inhaled 
vapors . 


Figure  5.  Effect 
of  increasing 
exposure  concen¬ 
tration  of  vapor- 
modifier  (abcissa) 
on  metabolic  rate 
of  the  vapor 
inhaled  at  constant 
concentration. 
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CONCLUSIONS 


We  have  presented  rationales  for  the  determination  of  metabolic 
constants  in  vivo  from  pulmonary  uptake  rate,  and  briefly  described 
the  procedures  and  some  applications. 

The  determination  of  the  extent  of  vapor  metabolism  from  the 
uptake  rate  has  an  advantage  over  making  this  determination  from 
excreted  metabolites,  in  that  the  effect  of  metabolite  distribution 
and  binding  in  the  body,  and  the  effect  of  renal  clearance,  are 
eliminated.  Since  air  sampling  is  a  noninvasive  procedure,  samples 
can  be  collected  continuously,  or  as  frequently  as  needed. 
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With  progressing  technology,  toxicologists  and  hygienists  confront  the  problem 
of  securing  safe  exposure  to  an  increasing  number  of  air  pollutants.  The  adverse 
biological  effect  of  pollutants,  like  the  therapeutic  effect  of  drugs,  is  related 
to  blood  concentration  and/or  time  Integral  of  concentration  in  the  target  organs. 
Passage  of  pollutants  from  the  environment  to  the  target  organ  has  the  same 
significance  as  migration  of  drugs  from  the  site  of  administration  to  the  target 
organ.  Pharmacokinetics  describing  the  transport  of  drugs  in  the  body  is  a  potent 
tool  for  designing  dosage  regimens  of  optimum  therapeutic  effect. ^  Methods  similar 
to  those  used  by  pharmacokineticists  can  be  employed  to  design  exposures  with  mini- 
mal  undesirable  biological  effects. 

Inhalation  administration  has  some  specific  characteristics :  Equilibration 

of  partial  pressures  of  Inhaled  vapor  in  the  body  and  in  ambient  air  is  the  driving 

(2  3) 

force  determining  uptake.  *  '  The  equilibration  rate  depends  on  pulmonary  ventila¬ 
tion,  tissue  perfusion,  and  on  solubility  and  clearance  cf  Inhaled  vapor.  The 
concentration  in  tissues  depends  on  exposure  concentration,  exposure  duration,  and 
equilibration  rate. 

Since  solubility  of  Inhaled  vapor  varies  for  different  tissues  according  to 
water  and  lipid  content,  and  since  cardiac  output  is  not  equally  distributed,  a 
multi -compartmental  model  is  needed  to  describe  uptake,  distribution,  and  elimina¬ 
tion  of  inhaled  vapors  (Figure  1). 


This  study  was  supported  by  AFOSR  762970  and  NIH  grant  ES01029-02. 
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In  the  model,  tissues  are  assigned  to  the  compartment  according  to  perfusion, 

ability  to  metabolize  the  inhaled  substance,  and  solubility  of  the  substance  in 
(4  5) 

the  tissue.  ’  'Lung  tissue,  functional  residual  air,  and  arterial  blood  form  the 
central  compartment  'LG' ,  in  which  pulmonary  uptake  and  clearance  take  place.  The 
partial  pressure  of  inhaled  vapor  equilibrates  with  four  peripheral  compartments. 
Well  perfused  tissues  form  two  peripheral  compartments:  BR-compartment  includes 
brain,  which  lacks  capability  to  metabolize  most  xenobiotics,  and  is  treated  as 
a  separate  compartment  because  of  its  biological  importance  and  the  toxic  effect 
of  many  vapors  and  gases  on  CNS.  VRG-compartment  includes  vessel  rich  tissues 
with  sites  of  vapor  metabolism  such  as  liver,  kidney,  glands,  heart,  and  tissues 
of  the  gastrointestinal  tract.  Less  perfused  tissues  are  also  pooled  in  two  peri¬ 
pheral  compartments,  according  to  lipid  content:  Muscles  and  skin  form  compartment 
'MG',  and  adipose  tissue  and  white  marrow  form  compartment  *FG'.  It  is  important 
to  treat  the  FG-compartment  separately,  since  dumping  of  lipid  soluble  vapors  in 
this  compartment  has  a  smoothing  effect  on  concentration  variation  in  other  tissues, 
caused  by  changes  in  exposure  concentration,  minute  ventilation,  and  exposure 

duration.  This  model  is  described  by  a  set  of  five  first-order  differential  equa¬ 
te) 

tions  linear  to  the  first  approximation.'  ' 

Mathematical  solution  of  this  model  is  available  if  the  model  is  pictured  as 

(4-9) 

an  electric  network  composed  of  conductances  and  capacitances:  '  In  figure  2 
'Z'  stands  for  exposure  concentration.  The  values  of  capacitances  ' C '  are  derived 
from  capacity  of  tissues  to  retain  the  vapor,  that  is ,  tissue  volumes  multiplied 
by  appropriate  tissue-air  partition  coefficients  'X'  at  37°C.  The  values  of  conduc¬ 
tances  SgR’^VRG  ’  SlG  an<*  ^FG  are  ^er1ve<^  ^rom  transportation  rates  of  vapor  from 
the  lung  to  the  tissues,  that  is,  blood  blow  times  blood-air  partition  coefficient 
(37°C).  Alveolar  ventilation  was  substituted  for  G^g.  'G^stands  for  clearance 
by  metabolism.  All  parameters  required  by  the  model  can  be  defined:  Physiological 
parameters  can  be  found  in  the  1  iterature^’^’^'^Partition  coefficients  can  be 


(2) 


easily  measured.  Metabolic  clearance  can  be  determined  from  pulmonary  uptake 
during  steady  state^0’^"*7^ 

The  determination  of  metabolic  rate  from  the  uptake  rate  has  an  advantage 
over  making  this  determination  from  excreted  metabolites,'  'in  that  the  effect 
of  metabolite  distribution  and  binding  in  the  body,  and  the  effect  of  renal 
clearance,  are  eliminated. 


After  sufficiently  long  exposure,  the  steady  state  is  reached  and  pulmonary 
uptake  rate  equals  clearance  rate.  If  clearance  does  not  take  place,  uptake  equals 
zero  and  ratios  of  tissue  concentrations  to  the  exposure  concentration  equal  the 
corresponding  partition  coefficients.  If  the  vapor  is  excreted  or  metabolized 
during  exposure,  the  ratio  of  concentrations  in  alveolar  air  and  in  tissue  to  expo¬ 
sure  concentration  is  smaller  than  corresponding  partition  coefficients.  The  de¬ 
viation  from  partition  coefficient  is  directly  related  to  clearance  and  indirectly 
related  to  the  flow  of  vapor  to  the  site  of  metabolism. 

The  uptake  rate  V  can  be  determined  from  difference  of  exposure  concentra¬ 


tion  *Cexp'  and  vapor  concentration  in  mixed  exhaled  air  ' Cgxh ' , 
'C  iv'  or  arterial  blood  * C  t * . 


alveolar  air 


u  =  (C  -  C  .  )  V 
'  exp  exh' 


(1) 


(2) 

(3) 


•  •  •  ,  • 
where  'V'  is  minute  ventilation,  'V  ^  '  is  alveolar  ventilation  (V^  *  2/3V)  and 

' x*  is  partition  coefficient. 


(3) 


Determination  of  uptake  rate  by  analysis  of  air  samples  (equations  1  and  2) 

has  the  advantage  that  sampling  of  mixed  exhaled  air  as  well  as  of  end  exhaled 

air  (alveolar  air)  can  be  done  frequently  without  imposing  stress  on  the  subject. 

M5  19-221 

However,  this  method  Is  limited  to  "cooperative  subjects",  such  as  men*  ’  ' 

or  animals  which  tolerate  a  face  maskf^*^  Sampling  of  arterial  blood  (equation 

(231 

3)  is  more  suitable  for  small  experimental  animals/  '  Anesthesia  or  any  drug 
administered  to  subjects  undergoing  the  exposure,  might  affect  the  metabolism  of 
inhaled  vapor. 

For  organic  solvents,  metabolism  is  the  main  excretory  pathway,  and  therefore 

«  • 

during  steady  state,  metabolic  rate  u^  =  u.  If  flow  rate  of  the  vapor  to  the 
site  of  metabolism  is  much  larger  than  metabolic  rate,  the  measured  clearance 
is  intrinsic  clearance  and 


(4) 


However,  for  most  vapors,  pulmonary  ventilation  and  tissue  perfusion  ( 1 F * )  affects 
the  metabolic  rate,  and  Gx  must  be  calculated: 


or 


or 


Cal  v  1 
Q  ^Xbl/air 


1  _  1  ^art  1\ 

G  "  A.  o  F 

x  bl/air 


_1 _ 

bl/air 


(5) 

(6) 

(7) 


(4) 


Gx  can  be  calculated  most  accurately  if  vapor  concentration  in  metabolizing  tissue 
during  steady  state  is  known: 


G 


x 


*ti s/alr 


(8) 


Perfusion  of  metabolic  rate  * F '  can  be  calculated  by  subtracting  equation  8  from 
equation  5,  6  or  7  (making  u  =  um) 

Gx  Is  a  constant  if  metabolism  follows  first  order  kinetics  (low  exposure 

concentration).  However,  metabolism,  like  all  enzymatic  reactions,  is  a  capacity 

limited  process.  Therefore  G  becomes  a  dependent  variable  of  tissue  concentration, 

(24) 

and  of  exposure  duration.  '  Gx  can  be  calculated  by  substituting  for  'u  1  in 
equation  8  from  Michael is-Menten  equation: 


tis/air 
Ct  1  s 


Vmax  Ctis  a 
Km+  Ctis 


Vmax  ^tis/air 
K  +  C' 


m 


'tis 


(9) 


K  and  V  can  be  determined  from  double  reciprocal  plot  of  uptake  versus  tissue 
m  max 

concentration  measured  during  steady  state  in  subjects  exposed  to  different  con¬ 
centrations  . 

The  non-linear  element  Gx  describing  metabolism  of  limited  capacity  further 
complicated  the  mathematical  solution  of  the  model.  However,  there  is  a  solution 
for  electric  network  with  non-linear  element  (25,26)  an{j  compUter  programs  are 
available. -j^g  pr0gran]S  calculate  time  course  of  voltages  in  capacitors 
and  currents  In  resistors.  The  voltages  'V'  are  directly  related  to  tissue  con¬ 
centrations  (C^  *  V  X^s/a-fr).  Currents  are  equivalent  to  uptake  rates  (or  to 
pulmonary  clearance  rate)  (in  GLQ),  metabolic  clearance  rates  (In  Gx), 
and  retention  rates  (in  GBR,  GVRQ,  GMQ,  GpQ),  Integrated  currents  represent 
amounts  retained,  metabolized  or  exhaled. 


(5) 


Examples  of  the  application  of  modeling  to  problems  of  industrial  toxicology 
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FIGURE  1 


FIVE  COMPARTMENTAL  MODEL  WITH  METABOLISM 
IN  VESSEL  RICH  COMPARTMENT 


FIGURE  3 


PREDICTED  CONCENTRATIONS  OF  BENZENE  IN  EXHALED  AIR  (CEXH)  OF  A  PERSON 
EXPOSED  TO  BENZENE  (8  HRS/DAY,  5  DAYS/WEEK)  FOR  6  WEEKS 


The  non-broken  line  counts  with  metabolic  clearance  *  3.6  1/min  (15,  27),  the 
broken  line  represents  the. hypothetical  situation  if  benzene  is  not  metabolized. 


Conclusion:  Metabolism  reduces  benzene  concentrations  in  exhaled  air.  The 
concentrations  at  the  end  of  the  week  are  larger  than  at  the  beginning  of  the 
week,  and  for  five  weeks,  rise  slightly;  on  the  sixth  week,  the  steady  state 
is  reached. 


(Reproduced  from  reference  4) 


FIGURE  4 


EFFECT  OF  SODY  8UILD  ON  EQUILIBRATION  OF  PARTIAL  PRESSURES  OF  INHALED 
VAPORS  IN  BRAIN  WITH  EXPOSURE  CONCENTRATION  DURING  8-HOUR  EXPOSURE 

NON-METABOLIZEO  METABOLIZED 


LOW  SOLUBLE  GAS 
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LIPID  SOLUBLE  VAPOR 
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The  partial  pressure  ratios  are  calculated  for  brain  of  a  normal  build  person 
(solid  lines),  a  slightly  obese  person  (dashed  lines)  and  a  slim  person  (dotted 
lines).  If  the  lines  coincide,  the  broad  solid  line  is  used, 


Conclusion:  Partial  pressure  equilibration  of  low  soluble  gas  in  brain  is  rapid, 
and  body  build  has  no  significant  effect  on  brain  concentration,  Equilibration 
of  lipid  soluble  vapors  Is  slow,  and  concentration  reached  In  brain  of  slim  person 
Is  much  higher  than  concentration  In  brain  of  obese  person.  If  the  Inhaled  vapor 
Is  metabolized,  the  brain  concentrations  are  reduced. 


(Reproduced  from  reference  5) 


FIGURE  5 


EFFECT  OF  FLUCTUATION  OF  EXPOSURE  CONCENTRATION  ON  CONCENTRATION  OF 

INHALED  VAPORS  IN  BRAIN 
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The  tissue  concentrations  are  presented  as  partial  pressure  ratio  of  tissue 
to  TVl.  Solid  lines  represent  8-hour  exposures  to  constant  concentration  (TLV). 
The  interrupted  lines  represent  examples  of  exposures  with  excursion  factors 
approved  by  ACGIH  (28):  1.25  (dashed  lines),  1.5  (dark  dotted  lines)  and  3 
(light  dotted  lines).  Excursions  for  each  hour  are  as  follows: 
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Conclusion:  Brain  concentrations  of  low  soluble  gas  (upper  graphs)  fluctuate 
with  exposure  concentration.  The  smoothing  effect  of  FG-compartment  on  brain 
concentrations  of  lipid  soluble  vapor  is  apparent  in  lower  graphs.  Metabolism 
diminishes  concentrations  reached  in  brain. 


(Reproduced  from  reference  5) 


FIGURE  6 


EFFECT  OF  STEL  ON  CONCENTRATIONS  OF  CYCLOPROPANE 
ANO  TOLUENE  IN  ARTERIAL  BLOOD 
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The  unbroken  lines  represent  blood  concentrations  during  8-hour  exposure  to 
constant  concentrations  of  cyclopropane  or  toluene.  The  broken  lines  represent 
blood  concentrations  at  the  end  of  a  single  15-minute  excursion,  with  two-fold 
or  three-fold  Increase  of  exposure  concentration.  The  values  are  plotted  for 
single  excursions  which  happen  any  time  during  exposure. 


Conclusion:  The  concentration  increase  in  alveolar  air  (and  in  tissues)  during 
short  excursion  is  always  smaller  than  the  excursion  factor  (increase  in  exposure 
concentration).  The  increase  depends  on  the  excursion  duration,  excursion  factor, 
duration  of  exposure  prior  to  excursion,  and  physical  and  chemical  properties  of 
inhaled  compound. 


r 1 bUKC  / 


EFFECT  OF  CAPACITY  LIMITED  METABOLISM  ON  TRICHLOROETHYLENE  CONCENTRATIONS 
IN  NON-METABOLIZING  ANO  METABOLIZING  TISSUES  DURING  STEADY  STATE 


NON-  MET ASOU2NG  METABOLIZING 


For  Rhesus  Monkey  ■  3.17±  0.73  mg/1  of  alveolar  air,  and  *  1.1-  0.12 
mg/min.  (16).  The  trichloroethylene  partial  pressure  ratios  tissue/exposure 
concentrations  are  plotted  against  exposure  concentrations.  The  curves  on  the 
left  refer  to  alveolar  air,  arterial  blood,  and  to  tissues  In  which  no  meta¬ 
bolism  of  trichloroethylene  occurs.  The  curves  on  the  right  refer  to  tissues 
and  venous  blood,  leaving  tissues  In  which  trichloroethylene  Is  metabolized. 
The  middle  line  Is  calculated  for  *  3.17  mg/1 ,  and  Vmax  »  1.1  mg/min.  The 
shaded  area  demonstrates  the  changes  caused  by  variation  of  Vmax  in  range  of 
2  S»0.  The  dotted  lines  demonstrate  the  changes  caused  by  variation  of  in 
range  of  2  S.D.  The  arrows  Indicate  the  direction  of  increasing  Vmax  or  Km 
respectively.  TLV  and  STEL  refer  to  threshold  limit  concentrations  (28),  MAC 
to  minimum  anesthetic  concentration. 


Conclusion:  The  concentrations  In  non-metabolizing  tissues  are  larger  than  in 
sites  of  metabolism.  At  low  exposure  concentrations,  the  partial  pressures  do 
not  equilibrate,  and  C^is  <Cexp  Xtls/alr*  When  exposure  concentration  approaches 
tissue  concentrations  start  to  rise  rapidly.  At  high  exposure  concentration, 
the  ratio  of  tissue  concentration  to  exposure  concentration  approaches  the  value 
of  corresponding  partition  coefficient.  The  tissue  concentrations  increase  with 
Increasing  and  decreasing  vmax- 


